RF dressed atoms beyond the linear Zeeman effect 
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We evaluate the impact that non-linear Zeeman shifts have on resonant radio-frequency dressed 
traps in an atom chip configuration. The degeneracy of the resonance between Zeeman levels is 
lifted at large intensities of a static field, modifying the spatial dependence of the atomic adiabatic 
potential. In this context we find effects which are important for the next generation of atom-chips 
with tight trapping. 

PACS numbers: 37.10.De, 37.10.Gh, 03.75.Be, 32.60.+i 



I. INTRODUCTION 

The use of radio-frequency fields (RF) for the manipu- 
lation of ultra-cold atomic samples is underpinning 
important developments in areas such as matter-wave in- 
tcrfcrometry [JJ ITHTT] which are beyond the original use 
of RF fields for evaporative cooling 12] . Nowadays, in 
combination with atom-chip technology |13j or optical 
lattices [TH [15] RF dressing is a well established tech- 
nique that allows us routinely to control and engineer 
atomic quantum states and potential landscapes on mi- 
cron scales [Ml [16] . Furthermore, RF dressing plays a 
central role in several proposals for extending the scope of 
functions and applications of ultra-cold atomic gases, in- 
cluding reduced dimensionality and connected geometries 
(ring and toroidal traps) [H IT7H22] , cooling and prob- 
ing of RF dressed atom traps [SI [531 121] > sub- wavelength 
tailoring of potentials [14] and transporting atoms in 
dressed atom traps [25] , 

Experimental realizations of RF dressed magnetic 
traps have worked in a range of static field intensities 
that produce linear Zeeman energy shifts [3j Q31 [23l I26T 
I5B] . Nevertheless, developments in near surface trap- 
ping/control and micro-fabrication |29j indicate that pro- 
duction of strong trapping configurations and sub-micron 
control will be soon on the agenda [301 SI] and a full de- 
scription of the atomic Zeeman shifts becomes relevant. 

One impressive application of RF dressing of magnetic 
traps is the miniaturized matter-wave interferometer for 
coherent spatial splitting and subsequent stable interfer- 
ence of matter waves on an atom chip [H [TT], [35]. In 
the interferometer, the potential landscape typically com- 
prises a double well in a transverse plane, accompanied 
by longitudinal weak trapping [JJ I35H54"] . It has been 
proposed that the double well potential can be helpful 
to study entanglement and squeezing phenomena, and to 
study phase coherence dynamics and many-body quan- 
tum physics [3H [35] . Here we establish the relevance of 
non-linear Zeeman shifts on the production of a double 
well potential in a typical atom chip configuration. 

For the purposes of this work, the potential landscape 
is produced by resonantly dressing a static magnetic field 
comprised of a quadrupole field with an offset field. Then 



the overall field components are 
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Here G is the magnetic field gradient near the quadrupole 
centre and Bo is the uniform offset field (see Fig. [TJ. In 
order to give quantitative results we focus on the ground 
state manifold of Rb87 as an example, and evaluate the 
dressed potentials for the \F = 2, mp = 2) state and also 
give some results for the \F — I, nip = —1) state. For 
simplicity, our analysis is restricted to small amplitudes 
of the RF field (B R p), such that beyond rotating wave 
approximation (RWA), effects can be ignored [33 j. This 
is also helpful to identify clearly the effects due to non- 
linearity of the energy shifts. To investigate the relevance 
in current experimental situations, we took parameters 
from recent experiments, i.e. we took G — 22.67 T/m, 
B = 1.0 or 3.0 G and B RF = 35.7 mG [34]: this enables 
us to quantify the effect that the non-linear Zeeman shifts 
have on the shape of dressed potential. 

In the following we first review the theory of the Zee- 
man effect in non-linear situations (section |n| and then 
apply the results to RF dressing in the strong field regime 
(section |III| . We report on the effects on both the vi- 
brational frequency of RF dressed atom traps and their 



locations, and then the paper concludes in section IV 



II. ZEEMAN SHIFT THEORY 

In the limit of slow atomic motion, the magnetic mo- 
ment of the atoms keeps its orientation relative to a spa- 
tially varying magnetic field, Bdc- The dynamics can 
be described by the Hamiltonian [351 133 



H = AI ■ J + n B \B D c\(gjJz - g'lh 



(2) 



where A is a measure of the hyperfine splitting, J and I 
are the spin and nuclear angular momentum operators, 
and the electronic and nuclear (/-factors are gj and g'j, 
respectively. The Bohr magneton is denoted by fis, as 
usual. 

As is well known, the operator F z = J z + I z commutes 
with the Hamiltonian ([2| and thus mp is a good quantum 
number (along with I and J). Since we consider Rb87 
with 7 = 3/2 and J = 1/2, the interaction mixes pairs of 
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Figure 1: (a) Schematic of a magnetic field quadrupole dis- 
tribution [Eq. |T])] and the RF coupling field. The horizontal 
double-headed arrow (red) indicates the polarization of the 
coupling RF field. The offset field Bo points out of the page. 

(b) Contours of a typical RF dressed double well potential. 

(c) Potential energy along y — (the x-axis) in panel (b). Pa- 
rameters characterizing the double well are defined as: well 
position x w , inter- well barrier height h and well frequency w. 



states \itif — 1/2, 1/2) and |mj? + 1/2, —1/2) in a mi, mj 
basis: \mj,mj). These two states have the same value 
of mp by construction. In this case, because J = 1/2, 
the diagonalization of Hamiltonian ^ reduces to 2 x 2 
matrix blocks and some uncoupled terms. This leads to 
a Breit-Rabi formula for the hyperfme energy spectrum 
of an alkali atom in a magnetic field [36 : 
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In the weak field regime, corresponding to a <C 1, the 
Zeeman shifts become linear in LB_dc| an d all these en- 
ergy levels [Eq. (J3j) and Eq. Q] are approximated by 
(see e.g. Ref. [37] ) 
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Here the signs in front of the square root contribution in 
Eq. ([3]) have been absorbed by the definition of the g- 
factor associated with the total angular momentum F — 
I+J. Because J = 1/2 there are only two possible values 
F = I + J = 2 and F = \I - J\ = 1: 
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Figure |2j[a) shows the atomic energy levels, correspond- 
ing to the upper hyperfme manifold, at different posi- 
tions in the field distribution Eq. |l]). In the case of 
RF dressed atom traps with Rb87 atoms in the trapping 
states mp = 1,2 follow such energy curves adiabatically 
and become trapped at positions of minimum field am- 
plitude pi]. 



III. DOUBLE WELL DRESSED POTENTIAL 

We consider a magnetic trap with static field distri- 
bution (JlJ, dressed by a uniform and linearly polarized 
oscillating magnetic field (RF) P] 2] . The coupling be- 
tween atomic states is dominated by the component of 
the RF field orthogonal to Bpjc [S], vanishing at positions 
where they are parallel (e.g. at positions of coordinates 
(x = 0,y)). A typical resulting double well potential is 
shown in Fig. [IJb). To characterize the potential energy, 
we evaluate the dressed energy along the x axis, where a 
double well appears as a consequence of the dressing. 

We define the z-axis parallel to direction of the static 
field, take the RF field polarized along the z-axis and 
apply the standard rotating wave approximation (RWA). 
Then, the Hamiltonian ([2| in a rotating frame becomes 



H 



1/2) 2 + 2m F a(gj + g]) + a 2 (gj + 9l ) 2 



AI-J + fi B \B DC \{gjJ z - g'lh) ± hw RF (j z + i z ) 
. Hb\Brf\ 
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(7) 



A(l + \ag\mp) 



where a = /ib\Bdc\/A- However, there are two uncou- 
pled states in the mj, raj basis: \I, 1/2) and | — I, —1/2). 



Because J = 1/2 the sign ± is chosen according to which 
(3) hyperfme manifold we are interested in: F = I + J or 
F = \I—J\. This is because each polarization component 
of the dressing field couples magnetic sublevels within 
a subspace with a given total angular momentum (see 
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Figure 2: (a) Zeeman shifts for the upper manifold (F = 2) of the ground state of Rb87, in a field distribution as Eq. 
Arrows indicate locations of resonance at low (left) and high (right) field intensities. Panels (b) and (c) show dressed potentials 
(solid) and bare states (dashed) for resonance located at (b) x w — 2 /j,m, with wrf/^-k — 2.31 MHz and (c) x w = 30 pirn with 
wrf/2tt = 14.4 MHz. Notice how, at large fields, several crossings of bare states appear. In all cases the field gradient is 
G = 22.67 T/m, the RF amplitude is B RF = 35.7 mG, and the offset field is B Q = 3 G. 



Appendix [A} . The diagonalization of Eq. ^ produces 
the dressed state energies displayed in Fig. [3](solid curve) 
for the magnetic field configuration of this paper. 

In the weak field regime and for RF frequencies much 
smaller than the hyperfine splitting (wrf A/H), the 
dressed energies for states of the upper hyperfine mani- 
fold are [1] 
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which produces, for the field distribution of Eq. |T]) and 
for weak-field seeker states [16], the dressed potentials 
shown in Fig. [3] (dashed line) . The double well is con- 
veniently described by three parameters: the well min- 
imum position (x w ), the height of the inter- well barrier 
(h) and the well frequency w which is defined through 
a harmonic approximation centred around the potential 
minimum [see Fig. file)]. To distinguish these proper- 
ties being evaluated according to Eq. ([8| or Eq. 
we denote the former quantities (i.e. from the weak field 
expression) with a superscript 0. 

The potential well's position and frequency are de- 
termined by the location of the resonant coupling, as 
suggested in Fig. 2lb)-(c), and by the intensity of the 
RF field. Figure [3 shows typical dressed potentials in a 
double-well regime (where we show only one of the two 
wells) . We compare both the dressed eigenvalue of Hamil- 
tonian Q , and Eq. ([8| for different RF frequencies wrf ■ 
A rather significant change in the well shape is seen for 
quite modest well separations for this field gradient. (The 



well separation is 2xq in Fig. [3j) 

With a fixed static field configuration, the distance be- 
tween the wells is controlled by the RF frequency wrf- 
At large distances, the non-linear character of Zeeman 
shifts impacts on the dressed potential, and Eq. ([sj) 
is no longer valid. In particular, the well frequency is 
significantly modified, as seen for mp — 2 in Fig. [4j 
For well separations of a few microns, the correction 
(Ax ° ) is about 5% (or 10%) for offset fields 

of 1 G (or 3 G), respectively. The relative effect on the 
frequency increases approximately linearly with the well 
separation, and can become a significant fraction of w 
for separations of tens of microns. The weakening of the 
well (or equivalently, the reduction of w), can be qualita- 
tively understood as a consequence of the appearance of 
multiple resonant positions due to non-linearity of Zee- 
man shifts [see Fig. [2jb)]: as the distance between bare 
state crossings increases, the dressed potential softens in 
comparison with the potential from a single resonant po- 
sition. This becomes one of the main effects of the non- 
linear Zeeman effect on dressed RF potentials. 

Figure [5] shows similar results for what is the F = 1 
manifold in the weak field regime. However, in this case 
we see that the trap is tightened by the non-linear Zeeman 
effect. The result is understood from the fact that gF, 
Eq. changes sign for F = 1. As a result the sequence 
of levels seen for F = 2 in Fig. [2^c) is reversed, which 
enables us to have tighter trapping as seen (inverted) 
at the bottom of the manifold of Fig. |2jc). Since the 
trapping on Rb87 F — 2 involves potentials twice as steep 
as those on Rb87 F = 1 we would not expect the trap 
to become tighter in F = 1 overall. (The factor two 
in steepness is because the two manifolds have the same 
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Figure 3: Comparison of non-linear and linear dressed potentials. There are two potential wells at ±xo and we show here 
just the well located at x — +xq with mp = 2: Red (solid) eigenvalue of Hamiltonian d7b. Blue (dashed) line is the adiabatic 
potential (jsj. The dressing frequencies and well positions are (c) 2.6 MHz, xo ~ 10 /jxa (b) 5.2 MHz, 30 /jxa and (a) 8.2 MHz, 
50 /xm. The quadrupole field gradient is G = 22.67 T/m [34] and the bias field is Bo = 3.0 G as in Fig. [2] In each of (a-c), 
wrf has been adjusted to produce a minimum at x^ — 1.0 /xm using the weak field expansion Eq. (10 1. 



magnitude of gp, which means a factor of two in the 
potentials because for F = 2 we can use mp = 2, while 
for F = 1 we can only use mp = — 1.) However, Figs. 
[4] and [5] show that, in absolute terms, F = 1 trap does 
become tighter than the F — 2 trap for quite modest 
values of the RF frequency. In essence, the weakening 
of the F = 2 trap by the non-linear Zeeman effect, and 
corresponding tightening of F — 1 , becomes sufficient for 
the trap frequency in F = 1 to become higher. 

The potential well's frequency for the dressed state 
\F = l,mp = —1) can be evaluated analytically using 
the solution of a three level system (see Appendix [B]) , 
and approximating the well's position by taking the po- 
sitions of the resonant couplings (large dots in Fig. [2jc)) 
and averaging them. The frequency obtained following 
this procedure coincides very closely with the numerical 
results shown in Fig. 5. 

In the case of F = 2, we find that, similarly, the 
well frequency corresponding to the dressed state \F — 
2,mp = 2) can be estimated by considering couplings be- 
tween the bare levels \F = 2, mp = —1), \F = 2, mp = 0) 
and \F = 2,mp — 1). This is shown in Fig. [4] (short- 
dash) and indicates that at high RF frequencies the cur- 
vature of the dressed energy is dominated by the cross- 
ings of just these three levels (see Fig. |2jc)). In the 
intermediate regime of dressing frequency, the coupling 
of all levels contribute significantly to the dressed state 
and the three level approximation is not valid. However, 
at low RF dressing frequencies, we observe that the lin- 
ear regime result coincides with the three level estimate 



of the well frequency scaled by v2 (dash and dash-dot 
line in Fig. [I| respectively). 

Figure [6] presents a comparison of the non-linear effects 
on the well position (Ax w = x w — x° ) and inter- well bar- 
rier height (Ah = h — h°). Again, these are calculated 
from the diagonalization of the Hamiltonian |7| and by 
comparing results to the weak field approximation, Eq. 
([8]). We see that the shift in the position of a potential 
well can be large in both the low RF frequency limit and 
the high RF frequency limit (see Fig. |6ja), inset). The 
shift at high RF frequency seems reasonable as the higher 
frequencies usually result in RF resonance in regions of 
stronger magnetic fields where the non-linear Zeeman ef- 
fect is expected to play a role. Note here that the vertical 
asymptotes in Fig. [6] correspond to the threshold were 
the RF frequency is just sufficient to excite transitions in 
the static field configuration, i.e. Hwrf — ¥ ^b9f\Bq\. 

To understand the shifts in position of the potential 
wells we look carefully at the issue of how we define the 
resonance location and the potential well location. Shifts 
of the well positions, Ax w — (x w — ie° ), can be under- 
stood by investigating where resonant RF coupling oc- 
curs, i.e, positions where Hwrp matches the energy sep- 
aration between levels with |Am^?| = 1. That is, 

Ellt ~ EZ m t±i = hw RF ( 9 ) 

which, in the weak field regime, reduces to a unique value 
of wrp given by 

Hb9f\B dc \ - tiw RF = 0. (10) 
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Figure 4: Well frequency as a function of RF frequency, obtained from the eigenvalue of Hamiltonian |7| (solid) and the 
linear approximation to the eigenvalues |8f (dashed), for (a) Bo = 1.0 G and (b) Bq = 3.0. Other parameters are as in Fig. [5] 
(including F — 2, = 2). The well frequency corresponding to a three level system and its value scaled by v2 are shown by 
the short-dash and dot-dash lines, respectively. 
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Figure 5: Well frequency as a function of RF frequency in the case F = 1, nif = —1. Results are obtained from the eigenvalue 
of Hamiltonian Q (solid) and the linear approximation to the eigenvalues (JsJ> (dashed), for (a) Bo = 1.0 G and (b) Bo = 3.0. 
Other parameters are as in Fig. [2] 



In the regime of weak magnetic field, multiple reso- 
nances occur at the same position [see Fig. [2^b)] , since 
the spacing between adjacent sub-levels is degenerate. In 
contrast, for large fields the degeneracy of transition fre- 
quencies is destroyed and the resonance locations become 
separated in space, as shown in Fig. [2jc). As in the treat- 
ment of the well frequencies above, a simple and good es- 
timate for the location of the potential well minimum x w 
can be found by taking an average of the positions where 
Eq. ([9| is satisfied for different (mj7,mf ± 1) pairs [dots 
in Fig. |2jc)]. This quantity, x w , has been evaluated us- 
ing both the procedure of Appendix B (Eq. |B7[ ), and 
by using the numerical solution of Eq. ^ itself. These 
results agree with the full solution of Eq. (JtJ) and arc 
presented as one curve in Fig. [6j 



IV. CONCLUSIONS 

Our calculations show that for a resonant RF dressed 
double well potential, the well frequency is a sensitive 
parameter to the non-linearity of Zeeman shifts. This 
is relevant for investigations of tunnelling processes in 
double well potentials, BEC interferometry, and phase 
coherence dynamics, since these phenomena are sensitive 
to the potential shape [3H[3S] . These modifications of the 
well frequency can be large enough to be tested exper- 
imentally via interferometry or atom cloud oscillations 

[am. 

We have also seen that in taking account of the non- 
linear Zeeman effect the sign of the g-factor is an im- 
portant consideration in the tightness of the resonantly 
RF dressed atom trap. In the case of Rb87 we found 
that for modest RF frequencies the F = 1 dressed trap 
could be tighter than the F = 2 dressed trap, which is 
quite against expectation in the linear regime. By using 
a model 3-level system, we calculate the well's frequency 
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Figure 6: (a) Main panel: Relative displacement of well po- 
sitions (Ax w /x° = (x w — x^ )/x° w ) versus Hwrf/ /isB , eval- 
uated via Eq. (M). On this scale, Ax w /x% are approximately 
independent of -Bo- Inset: Well position displacement for off- 
set field Bo = I G (dashed) and Bo = 3 G (short-dashed) (b) 
Modification of the double well barrier height (Ah — h — h ) 
for Bo = I G (solid) and Bo = 3 G (dashed) respectively, as 
function of the RF frequency. In all cases Brf = 35.7 mG 
and G = 22.67 T/m [34]. Dashed vertical lines indicate the 
limiting frequencies at fiBgFBo/h. 



for F = 1. In case of the F — 2 manifold, this model 
works well in a regime of non-linear Zeeman shifts. 

Finally, we note that near the limiting frequency set 
by the static offset field (see Fig. |6|, the double well sep- 
aration and frequency of the dressed potential are quite 
sensitive to the dressing frequency w RF (see also Figs. [4] 
and [5]). This regime is relevant for double wells with sub- 
micron separations, a situation likely to occur in near sur- 
face trapping configurations |80j . Our work shows that to 
achieve stable configurations with small well separations, 
good control of the RF frequency is needed. However, the 
main result of this work is that the vibrational frequency 
can be sensitive to the chosen RF frequency because of 
the non-linear nature of Zeeman shifts. 
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Appendix A: RF dressing in the strong field regime 

An alkali atom interacting with a static magnetic field 
plus an orthogonal RF dressing field is described by the 
Hamiltonian: 

H = AI ■ J + v B \B D c\(gjJ z - g'lh) 

+HB\B RF \(gjJ x - g'jl x ) cos(w RF t) (AI) 

where the first term arises from the hyperfine interac- 
tion. Assuming that the hyperfine coupling is stronger 
than the interaction with the static field, the state space 
can be split into a direct sum of spaces with two angular 
momentum operators (which has spin I — J) and F^ 
(which has spin I + J). Left and right circular polar- 
ization components of a linearly polarized dressing field 
couple magnetic levels within only one of the hyperfine 
subspaces. Thus, the relevant component in each case is 
selected by an observer rotating in an appropriate sense. 
With this in mind, the unitary transformation between 
lab and rotating frames can be defined by [35]: 



U - exp(-iw RF t(F^ - Fj)) , 



(A2) 



where FJ and i*j work in subspaces of F as described 
above. Then, in the rotating frame, the time evolution 
follows the Schrodinger Eq. ihdt\ip) — H'\ip) with 



H 1 = U^HU - iMJ^U. 



(A3) 



The hyperfine term and the interaction with the static 
field are invariant under the unitary transformation Eq. 
( A2 1 . The interaction with the dressing field transforms 



into: 



Hi 



RF 



HbB R f,x 
2 



(cos2w RF t x =p sh^u;^ y) ■ (gjj - g'jl) 

-{gjJx-g'M- (A4) 



Neglecting the rapidly rotating field (rotating with an- 
gular velocity of 2w RF ), the total Hamiltonian becomes 
the expression given in Eq. ([7]). 



Appendix B: Frequency and position of strong-field 
potential wells 

Consider the dressed energy as function of the static 
magnetic field, as given, e.g. in Eq. and produce a 
Taylor expansion around the minimum occurring at the 
magnetic field \B DC , m in\, 



E(\B DC \) = E + ^(\B DC \ 



\BDC.min\) + 



with /?, 



d 2 E(\B DC \ 



d\B 



-DC I 



(Bl) 



(B2) 



\B DC ,mi n \ 
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The well's frequency is parametrized by a quadratic de- 
pendence of the dressed energy with the distance to the 
position of minimum 



where 



. . „ raw , 

E{x)=E + ^—{x 



\BDC,min\ — \Bdc(%w)\- 



(B3) 



(B4) 



In the case of a quadrupole magnetic field of gradi 
ent G plus offset field, where the field is \Brjc( x )\ = 



Bq , the well's frequency in terms of expansion 



Eq. (Bl) is given by: 



mA 



Bo 



\B 



DC. rain \ 



(B5) 



Then in the regime of linear Zecman shifts, and after a 
Taylor expansion to second order of Eq. {[sj) we obtain 



\B 



hw 



RF 



DC.min \ 



P = 



SAg F 



RF I 



(B6) 



For this to be valid, the minimum should occur at a mag- 
netic field such that /j,b\Bdc. //</// I < 2A. 

For more intense fields, where the multi-level crossing 
degeneracy is lifted, but levels can still grouped in man- 
ifolds F = I + J and F = \I — J\, \Brjc,min\ can be 
obtained as an average of the crossing points between 
consecutive mp levels. The field at which such crossings 
occur can be obtained analytically by solving Eq. ([9]), 
expanding the Zeeman shifted energies (|3| up to second 
order in [IbB/A. This procedure gives us 



with 



G F = T- 



9.1 - 9i 



1+1/2 V2(7 + 1/2) 



(B8) 



In Eqs. (B7)-(B8) the upper and lower signs are chosen 



according to F = I + J and F = \I — J\ respectively. 

Analytic expressions for the dressed energy can only 
be obtained in simple cases. For our example of Rb87 
in its ground state, after ignoring couplings between the 
F = I + J and F = \I — J\, the dressed energies of the 
F = 1 manifold are: 



Ei(\B DC \) = -- + 3 cos 



with 



+A (-I + 2£|< / + »' 5 >) < B<J » 

0, 2-7T, Att corresponding to i — 1, 2, 3, and 
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D 



= -(Et 



z,- 



E, 



z,- 



E 



Z,-\ 



Ef'~E z 'i + E^'~Eq' 



E Z _\E Z ^ 



2d 2 + Hwrf(E Z ' — EZ\\ — Hwuf) 

E Z J-) - E z - E z >- E Z J- 
z, 



z,- 



E = 
R = 

Q = 



d 2 (E z 



-hw RF E z ' (E z ' 



ET 



hw RF ) 



(9CD-27E-2C 3 )/54: 

{3D-C 2 )/9 

arccos ( R / \f— Q 3 ) , 



(BIO) 



where d = MBg f FgF (m F ±l\J x +I x \m F ) and the Zeeman 
shifted levels are given by Eq. 



m F =F-l 



\B 



DC.min \ 



2F lH 
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9F 



9f "F 4(1 — 2m F )G F 



flW f. 



2(1 - 2m F )G F 
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